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Abstract 

We study the instability of a superfluid flow through a constriction in three spatial dimensions. 
We consider a Bose-Einstein condensate at zero temperature in two different geometries: a straight 
waveguide and a torus. The constriction consists of a broad, repulsive penetrable barrier. In the 
hydrodynamic regime, we find that the flow becomes unstable as soon as the velocity at the 
classical (Thomas-Fermi) surface equals the sound speed inside the constriction. At this critical 
point, vortex rings enter inside the bulk region of the cloud. The nucleation and dynamics scenario 
is strongly affected by the presence of asymmetries in the velocity and density of the background 
condensate flow. 
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Introduction. 

Even when flowing past an obstacle, a superfluid can be stationary, since its excitations 
appear only under certain conditions For instance, a homogeneous Bose-Einstein con- 
densate (BEC) of density n running through a weak obstacle, whose presence only slightly 
changes the cloud chemical potential, is energetically unstable towards phononic excitations 
only above a critical velocity v c equal to the sound speed c oc y/n, according to the Landau 
criterion. On the other hand, the critical velocity and the nature of the excitations change in 
the presence of a strong obstacle that considerably perturbs the homogeneous flow. Helium 
flowing through narrow channels becomes unstable at velocities much lower than predicted 
by the Landau criterion, which this time would correspond to rotonic excitations, since, as 
first predicted by Feynman, vortices are nucleated inside the channel and carry energy away 
from the superfluid through the phase-slip mechanism 2|-[5j. Feynman's energetics argument 
is analogous to the Landau criterion, and is based on the choice of vortex rings, instead of 
rotons, as the excitations coming at the critical velocity. His conjecture has been refined by 
Varoquaux [5], who used half-rings as the excitations to be nucleated inside the channel. In 
general, the type of excitation and the critical velocity depend on the flow configuration. 

In contrast to nearly-incompressible fluids like helium, atomic BECs, due to an exquisite 
control over the confining potential, offer a wide choice of flow geometries and obstacles, and 
also permit the modification of the effective dimensionality. Imaging techniques also provide 
a means of addressing the system on the healing length scale, at which vortex dynamics takes 
place. The superfluid critical velocity in an elongated BEC can be studied for instance by 
sweeping a laser beam through the cloud jfj. At the critical velocity, both vortices in an 
elongated trap [7f and solitons for tight confinement transverse to the flow [8j have been 
experimentally observed. Another way to study the stability of the superflow would be to 
perturb a persistent current, which has been experimentally demonstrated with a BEC in a 
torus [9(1 . Apart from being an hallmark of superfluidity, BECs in toroidal geometries and the 
study of their critical velocities have very important technological applications for sensing 
devices, like SQUIDs already realized with superfluid helium and superconductors. BECs 
are often very well described by the mean- field Gross-Pitaevskii equation (GPE), which 
provides, unlike for helium, a reliable theoretical model to study the instability mechanism 
and its dynamics. The nucleation of a vortex-antivortex pair at the critical velocity has 
been studied numerically in two and three dimensions with an obstacle moving through a 
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uniform condensate [lo|, and in a two-dimensional waveguide with an orifice [llj]. Flow 
dissipation caused by the formation of solitons in one dimension has been investigated in 
12]. In a previous work [l3], we studied the scenario of superfluid flow dissipation in an 



effective two-dimensional toroidal geometry. 

In this manuscript, we study the critical velocity and superfluid flow dissipation mecha- 
nism in a three-dimensional constricted flow configuration, where the size of the cloud along 



all the three directions is much larger than the healing length £ = h/ y/2mgn. We consider 
a subsonic flow of a zero-temperature BEC in two different geometries: a wave guide with 
periodic boundary conditions, which can mimic an elongated cloud along the flow direction, 
and a torus. The unstable regime is reached by raising a repulsive penetrable barrier per- 
pendicular to the flow. This barrier is broader than the cloud dimension and extends over 
a few (typically 5 to 10) healing lengths along the flow direction. In this way, we create a 
constriction for the flow in the barrier region, as shown in Fig. [TJ Starting from a stationary 
flow with a given velocity, the barrier is adiabatically raised during the dynamical evolution 
until the instability sets in, and we observe vortex rings penetrating the cloud. In this way, 
the superflow energy is dissipated through the phase-slip mechanism j^j: vortices are able 
to take energy away from the flow since the velocity drops by a quantized amount in a 
region crossed by a vortex core 3]. The two geometries under study present significative 
differences in the vortex ring nucleation and dynamics. We find that vortex rings, which can 
find a stationary configuration after entering an axially symmetric waveguide, are instead 
always transient in the torus, and, more generally, as soon as the axial symmetry in the 
direction of flow is broken. 

Formulation and Computation 

In order to evolve these various cases, we solve a GPE ([T]) for three spatial dimensions in 
scaled form as 

dif)[r,t) 

where length, time, and energy are given in units of d Q = [^-] 5 , l/w G , and ftw respectively 
for a representative harmonic frequency u Q that characterizes the trap. V(r, t) is the external 
potential and g = with a the inter-particle s-wave scattering length and m the atomic 
mass. if)(r,t) is the condensate wavefunction normalized to the total number of particles 
N. The external potential has components associated with the trapping and the barrier 



±V 2 + V(r,t)+g\ij\ 2 



i>(r,t) (1) 
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FIG. 1: (color online). Constriction configuration. The light-gray surface corresponds to the 
classical (Thomas-Fermi) surface of the cloud. The dark-gray surface shows an isosurface of the 
barrier potential used to create the constriction. 

potentials of the form 

V(r,t) = V tr (r) + V b (r,t). (2) 

In both cases, we find the ground state of the system with Vb = 0. 

The waveguide geometry is implemented by choosing the trapping potential 

V tT (r) = 1 -[x 2 + 1 2 z 2 }, (3) 

with uj = uj x = 30 x 27rHz, 7 = u z /u x , and periodic boundary conditions along the flow 
direction y. We considered three different values of the 7, namely 1, 1.05, and 1.2, which 
correspond, respectively, to a ground state chemical potential \i = 11.7, 12.0, 16.5 with 
N = 3 x 10 5 87 Rb atoms and a nonlinear scaling value g N = 10134 . 

The toroidal geometry is implemented by choosing the trapping potential 

V tT (r) = \[a 2 x 2 + (3 2 y 2 + z 2 ] + V c e~ 2 ^ , (4) 

where a = uj x /u z and j3 = uj y /u z . We take a — j3 — 0.5 with u Q = u z = 25 x 2nB.z and 
form the torus by including a core potential with parameters V c = 144, and a c = 1.88 with 
p 2 = x 2 + y 2 . The ground state chemical potential is fi — 7.6 for iV = 2.5aT0 5 23 Na atoms, 
corresponding to a nonlinear scaling value g N = 2028. 

During the dynamical evolution, we use a time-dependent barrier potential of the form: 

V b (r, t) = f(t)V s VUx)V hy (y)V hz (z)/8, (5) 
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with f(t) = t/t z (f(t) = 1 for t > t T ) and V bx = tanh( *-^+*° ) + tanh ^+g-^ ). Here it, 
is the x-shift of the center of the barrier while its width is w x ~ 2xq. Vb y (y) and Vb z (z) have 
the same form as Vt, x . The final height of the barrier is as long as xq, yo, z 3> b s . 

Before starting the dynamics, we put the condensate in motion by imprinting an ap- 
propriate spatially dependent phase 8(r) on the wavefunction. In the waveguide case, 
9(r) = mvy/h generates a uniform flow of velocity v along the y direction. In the torus, 
6{r) = ml(j)/h, with / integer, generates a tangential flow of speed v(r) = l/p, where 
= arctan(y/x). 

The results presented in this manuscript are obtained by solving numerically the GPE. 
After finding the ground state with an imaginary-time propagation, we perform the dy- 
namical evolution in real time. For all the simulations performed we used a finite-element 
discrete variable representation (DVR)[22]. In the waveguide geometry, we employed a split- 
operator method, with a Fast-Fourier- Transform algorithm used for the kinetic part of the 
evolution, while for the simulations in the torus we used a Real Space Product Formula 



221 ] . The spatial grid for the waveguide calculations was n x = 60, n y = 120, n z = 60 (or 



n x = 180, n y = 120, n z = 180 for Figs. Hand EJ) with box sizes L x = 12, L y = 24, L z = 12 
for the axially symmetric case (see below), n x = 90, n y = 120, n z = 80 with box sizes 
L x = 18, L y = 24, L z = 16 for the axially asymmetric case, while the time-step was 
dt = 10~ 4 . In the torus simulations, we used a spatial grid consisting of 80 elements in 
each dimension with 4 DVR Gauss- Legendre functions in each element spanning cubic box 
lengths of [—12, 12], [—12, 12], and [—10, 10] in the x, y, and ^-directions respectively. This 
choice gave 241 grid points in each direction. We also tested the convergence with 5 basis 
functions and 321 points with only a few percent change in basic quantities such as energies, 
momentum, and positions. 
Criterion for instability 

In the hydrodynamic regime, when the healing length £ is much smaller than all other 
length scales (the smallest of which is the barrier width along the flow direction), we find 
that the onset of instability coincides with a simple condition: inside the barrier region, the 
local fluid velocity at the classical (Thomas-Fermi) surface of the cloud equals the sound 
speed, as depicted in the inset of Fig. [2 It is important to precisely define the sound speed 
which sets the threshold for the critical velocity. The latter is the Bogoliubov sound speed, 
calculated inside the barrier region, for the low-lying modes propagating along the flow, 
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FIG. 2: (color online). Ratio of the higher local fluid velocity at the Thomas-Fermi radius, v(Rtf), 
to the sound speed c inside the barrier, as a function of the barrier height. The red-shaded area 
correspond to the critical point plus uncertainty. The results are obtained for the cylindrically 
symmetric waveguide 7 = 1 and initial flow velocity v = 1.05. The inset shows a sketch of the 
behavior of the local fluid velocity along a radial cut inside the barrier region. The red solid line 
indicates the value of the sound speed c, which with transverse harmonic confinement is just the 
average of the local sound speed, c(r±) (red dashed line) over the transverse plane. The black(blue) 
solid line corresponds to a subcritical(critical) condition. 

taken as if the system was homogeneous in this direction. In the case of a waveguide with a 
harmonic transverse confinement, and within the Thomas- Fermi approximation, this sound 
speed is simply the average of the local sound speed on a plane perpendicular to the flow 
c = c(0) / \/2, where c(0) is the local sound speed at the center of the transverse harmonic 
trap Q. 

We verified this numerically in the waveguide shown in Fig. [2j In a waveguide 

with cylindrical symmetry, this happens simultaneously at points on a circle perpendicular 
to the flow, from where a vortex ring will then enter. In the torus, due to a higher flow 
speed at the inner edge of the annulus, the critical condition is reached first on the interior 
of the cloud, with the consequence that vortex rings, if ever formed, must be asymmetric, 
as we discuss below. 
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Slightly below the critical point, we observe vortices getting closer to the edges of the 
cloud but failing to enter. The presence of these "ghost" vortices [13I, Q suggests that a 
pre- instability is triggered at the edges of the condensate, but, since the vortices do not 
enter the bulk region, it is not sufficient to dissipate the superflow. 

In a previous work we verified the same condition for instability to hold for a two- 
dimensional toroidal BEC. The results reported here represent an extension to a three- 
dimensional case. As discussed above, this criterion allows for a simple understanding of 
the details of vortex penetration dynamics, based on the observation that vortex cores enter 
first where the critical condition if first reached. 

Instability dynamics in the waveguide 

In analogy with the standard situation where a superfluid flows through a channel, we 
consider a waveguide with periodic boundary conditions along the flow direction y, and the 
barrier creating the equivalent of the channel. This situation can be realized experimentally 
with a BEC inside an elongated trap along the flow direction. As stated above, before raising 
the barrier in the simulations, a stationary flow is created by imprinting a phase mvy/h on 
the condensate wavefunction, where v is the constant flow speed. This corresponds in the 
experiment to sweeping the barrier across the cloud at a constant velocity. We consider a 
symmetric harmonic confinement in the transverse x — z plane, such that d x = ^Jh/muj x and 



d z = ^h/mu z are both sufficiently larger than the bulk healing length, so that the effect 
of transverse degrees of freedom comes into play, giving rise to a fully three-dimensional 



dynamics. Such a setup has been used in [16J, where, moreover, the possibility of creating 



a penetrable repulsive barrier with a control over a length comparable to the healing length 
has been demonstrated. 

In this configuration, as soon as the critical barrier height is reached, a vortex ring 
detaches from the system boundaries and starts shrinking into the cloud inside the barrier 
region, as shown in Fig. [31 a). In the figure, black points indicate the position of vortex cores 
while the gray surface corresponds to the Thomas- Fermi surface of the cloud. Detection of 
vortex cores is performed using a plaquette method, described in [21]. 

Depending on the initial flow velocity, waveguide transverse section, and barrier height 
at the critical point, the ring attains a certain radius and velocity with which it propagates 
in the flow direction in a stable fashion, as long as axial symmetry is preserved, as depicted 
in Fig. [31 c) and d). Here the vortex ring eventually propagates at the speed u r = 0.65, and 
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FIG. 3: (color online). Four subsequent stages of the vortex ring penetration in the waveguide. 
The gray surface indicates the position of the Thomas- Fermi surface of the condensate. Black dots 
show the position of the vortex cores. Here the waveguide is axially symmetric 7 = 1, the initial 
flow velocity is v = 1.05, and the final barrier height is V s = 0.17 '/i. a) The ring is shrinking around 
the cloud in barrier region, still outside the Thomas-Fermi borders of the cloud, b) The ring has 
just entered the cloud, c) The ring has shrunk to its final size and is already outside the barrier 
region, moving along the flow direction, d) The ring has moved far from the constriction region, 
with a constant speed and radius. 

a radius R = 2.9. 

For sufficiently strong barriers, the ring shrinks to a point crossing the transverse section, 
thereby annihilating and completing a full single phase-slip, as shown in Fig. HI After this 
process, the velocity has dropped everywhere by the same quantized amount. We observe 
that such an event is the three-dimensional analogue of vortex-anti-vortex annihilation in 
two dimensions fl3 |. It is interesting to analyze the ring annihilation event in more detail. 
In Fig. |5] a) and b), respectively, we show the position of the vortex cores toghether with 
the density on a plane parallel to the flow direction, at a time just before and just after the 
ring has shrunk to a point. In Fig. E]c), the density on the same plane is plotted at four 
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FIG. 4: (color online). Four subsequent stages of the vortex ring annihilation in the waveguide. 
The gray surface indicates the position of the Thomas-Fermi surface of the condensate. Black 
dots show the position of the vortex cores. Here the waveguide is axially symmetric 7 = 1, the 
initial flow velocity is v = 0.52, and the final barrier height is V s = 0.94/i. a) The ring is shrinking 
around the cloud in the barrier region, still outside the classical borders of the cloud, b) The ring 
has entered the cloud, c) The ring is about to shrink completely and annihilate, d) The ring has 
annihilated and no vortex core is now inside the Thomas-Fermi surface. 

subsequent times after the annihilation. The details of self-annihilation process which we 
observe are consistent with the previous studies of axisymmetric vortex ring solutions 



23- 



271 ] . Namely, if we consider the position of the vortex cores, we see that the points of phase 
singularity form a loop which shrinks inside the constriction, Fig. [5] a), and whose radius 
eventually becomes zero, Fig. [5]b). At this moment, we see the zero-density core being filled 
with atoms, thereby transforming into a density depression, which further propagates out of 
the constriction as a rarefaction pulse, as shown in Fig. c). Finally, this rarefaction pulse 
decays into sound. 

Both in Figs. [3] and HJ we see that, at least initially, the vortex cores are subjected to 
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t = 4.98 




t = 5.40 




FIG. 5: (color online). Details of the ring self-annihilation event studied in Fig.U at six subsequent 
times. In a) and b), left panel, the gray surface indicates the position of the Thomas-Fermi surface 
of the condensate, while black dots show the position of the vortex cores. In a) and b), right 
panel, and c), the density on a z = plane parallel to the flow direction is shown. A very small 
loop structure of vortex cores, in a), shrinks to a point and has disappeared in b). The ring has 
transformed into a rarefaction pulse, whose propagation and decay into sound appears in c). 



an essentially radial motion, giving rise to the shrinking of the vortex ring. This behavior 
can be understood by considering what contributes to the motion of vortex cores in non- 
homogeneous systems 191 ] . Let us consider a single vortex core located at some position 



x c . Its velocity would be the sum of two terms: i) the background flow velocity at x c , 
and ii) a term perpendicular to the gradient of the density at x c . Both contributions must 
be calculated as if the vortex was not present. Since the relative weight of term ii) is 
proportional to the value of the healing length calculated at x c , in high density regions, a 
vortex core will move mainly with the background superfluid velocity, while in low density 
regions, it will move mainly due to the gradient of the density. Therefore, when the vortex 
ring is inside the constriction, where the density is low, and especially when it is close to 
the Thomas- Fermi surface, it will mainly move perpendicular to the gradient of the density, 
whose main contribution comes from the density modulation induced by the barrier along 
the flow direction. This results in an essentially radial motion of the cores, and thus into 
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FIG. 6: (color online). Three subsequent stages of the vortex ring annihilation in the torus. The 
gray surface indicates the position of the Thomas-Fermi surface of the condensate. Black dots show 
the position of the vortex cores. Here the initial circulation is I = 1, and the final barrier height 
is V s = 0.5//. a) The ring is shrinking around the cloud in barrier region, still outside the classical 
borders of the cloud, b) The ring has entered the cloud, c) The ring is about to shrink completely 
and annihilate. 

the shrinking of the ring. 

Instability dynamics in the torus 

In the toroidal geometry, the scenario is richer since the tangential flow velocity at a 
given total angular momentum, due to the conservation of quantized circulation, decreases 
like 1/r, where r is the distance from the center of the torus. This introduces an asymmetry 
between the inner and the outer edges of the cloud, which is not present in the waveguide 
case. Consequences of this asymmetry on vortex nucleation in two dimensions are discussed 
in [3]. In three-dimensional configurations, the result is that vortex rings are transient 
features in a toroidal geometry. 

The ring vortex either breaks or shrinks to a point and annihilates for sufficiently strong 
barriers. While the annihilation process, shown in Fig. |6j is very similar to the one taking 
place in a waveguide, the vortex ring breaking mechanism reflects instead more evidently 
the asymmetry in the velocity field. 
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FIG. 7: (color online). Four subsequent stages of the vortex ring breaking in the torus. The gray 
surface indicates the position of the Thomas-Fermi surface of the condensate. Black dots show 
the position of the vortex cores. Insets show the top and side views. Here the initial circulation 
is I = 4, and the final barrier height is V s = 0.2/i. a) The vortex ring is bending to form a right 
angle, b) The vortex ring has formed a right angle whose vertex is close to a vortex line coming 
frome the center of the torus, c) The vortex ring and line have just reconnected: a vortex line and 
a portion of a ring vortex are now inside the Thomas-Fermi surface, d) The vortex line and the 
ring have moved apart. 

As shown in Fig. [7J the ring is strongly deformed since the inner part moves faster (the 
velocity of a vortex core, when the healing length is much smaller than the length scale of 
density variation, is essentially given by the background flow velocity (l9j). The deformation 
increases in time up to the point at which the ring bends in on itself, forming a right angle 
at whose vertex a kink is present, as shown in Fig. b). Meanwhile, a vortex line coming 
from the inner core of the torus approaches the vertex, also forming a kink in correpondence 
to the latter. Eventually, the line and the ring connect, joining each other at the position 
of the kinks. Such an event produces a vortex line plus a new ring vortex, appearing in 
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FIG. 8: (color online). Three subsequent stages of the vortex ring formation in the torus. The 
gray surface indicates the position of the Thomas-Fermi surface of the condensate. Black dots show 
the position of the vortex cores. Insets show the side view. Here the initial circulation is I = 4, 
and the final barrier height is V s = 0.2//. a) The vortex line is bending around the cloud in barrier 
region, partially outside the Thomas-Fermi borders of the cloud, b) The vortex line has developed 
two kinks, c) The vortex ring has just formed. 

Fig. c). Thus, the breaking is actually a vortex reconnection. Vortex reconnections play 
an important role in turbulent scenarios, being a very efficient mechanism for increasing the 



number of vortices in the system 20j . When the barrier is not strong enough to make both 
the edges of the annulus unstable, we also observed cases in which the ring is not formed at 
all, and a strongly bent vortex line enters the cloud from its inner edge, to circulate around 
the torus. The fact that the vortex line enters the inner edge of the annulus is due to the 
above mentioned asymmetry in the velocity field, decaying as the inverse of the distance 
from the center of the torus, which makes the instability set in there first, according to the 
criterion discussed in section „ 

The formation of a vortex ring in the torus is very interesting, since it can be seen as a 
reconnection of a vortex line with itself. As shown in Fig. [HJ a bended vortex line is always 
present at first. While the bending increases, the vortex line develops two sharp kinks whose 
tips get closer to each other, up to when they join, thereby cutting the original line into a 
vortex ring plus a yet another line. The latter is then reabsorbed at the system's boundary. 

Connection between effective two- and three-dimensional geometries 
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FIG. 9: (color online). Four subsequent stages of the vortex ring dynamics in the axially asym- 
metric waveguide. The gray surface indicates the position of the Thomas-Fermi surface of the 
condensate. Black dots show the position of the vortex cores. Insets show the front and side view. 
Here the waveguide is non-axially symmetric 7 = 1.2, the initial flow velocity is v = 1.05, and the 
final barrier height is V s = 0.13/i. a) A partially-ghost ring vortex has formed, b) The ring vortex 
is strongly deformed, c) The ring vortex has broken up leaving two vortex lines, d) The vortex 
lines have joined back to form a new ring vortex. 



When the condensate is effectively in two dimensions, meaning that the cloud size along 
the third direction is comparable with the healing lentgth £, vortex lines oriented along the 
direction of tighter confinement can become the preferred excitations with respect to vortex 



rings [18|J. In studying a condensate inside an effective two-dimensional toridal trap (the sys- 
tem's size along the axis of the torus was about one healing length) [13|, we indeed observed 
that superfluid flow dissipation took place through the formation of vortex lines entering the 
cloud. For the trap configurations considered in this manuscript, we have observed vortex 
ring formation in the very low density regions outside the condensate classical Thomas- Fermi 
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surface. As suggested above, we could call these ghost [14|| vortices since they are not visible 
from the condensate density profile. In some cases, as discussed above, these ghost vortex 
rings are able to enter completely the classical surface of the cloud, transforming into what 
we should then call "real" ring vortices. 

In the crossover between effective two- and three-dimensional regimes, moving from a 
scenario in which only vortex lines are present to one in which real vortex rings come into 
play, it is reasonable to expect an intermediate regime in which ghost vortex rings are formed, 
but the condensate is sufficiently squashed along the third dimension that the full vortex 
loop is not able to enter the cloud's classical surface. In this regime, superfluid would be 
dissipated by vortex rings which are partly real and partly ghost, appearing as simple bent 
vortex lines in the density profile. 

An example of such situation is given in Fig. [HI where the waveguide is non-axially 
symmetric about the flow direction (7 = 1.2). After the instability sets in, a ghost vortex 
ring forms and shrinks around the cloud in the barrier region up to when we observe a full- 
fledged ring vortex, which is only partially inside the Thomas- Fermi surface of the condensate 
(Fig. |9^)). However, since the part of the ring which is inside the Thomas- Fermi surface 
moves with a larger speed along the flow direction with respect to the part which remains 
outside, the ring vortex is soon deformed (Fig. IHb)), and eventually breaks up (Fig. IHfc))- 
The vortex cores located inside the Thomas-Fermi surface of the cloud move, to a good 
approximation, along with the background velocity field. On the other hand, cores in the 
low density region outside the surface of the cloud essentially do not feel the background 
velocity and move along the flow direction only because of the presence of transverse density 
gradients. After the vortex ring breaks up, the two lines move downstream and continue to 
deform, and eventually re-join to form a vortex ring (Fig. MX))- The latter undergoes the 
same deformation described above, leading to another break up. 

We also verified that, even in a very slightly non-axially symmetric waveguide, ring vor- 
tices eventually break up. With 7 = 1.05, the deformation is created more slowly with 
respect to the 7 = 1.2 case, but the ring breaks up anyway, though at a later time. 

Conclusions 

We studied the superfluid instability condition and dynamics in the presence of a three- 
dimensional constriction across the flow. We solved the GPE, which would accurately de- 
scribe an ultracold dilute Bose gas, but would anyway provide a useful insight for the under- 



15 



standing of the instability in strongly correlated superfruids like Helium II. We considered 
two different setups: a waveguide, which could be experimentally realized by trapping an 
ultracold Bose gas inside a cigar-shaped trap, and a torus. We found that the instability 
sets in when the fluid velocity at the Thomas-Fermi surface of the cloud equals the sound 
speed inside the constriction region. This critical condition explains the strong dependence 
of the instability dynamics on the geometry considered, and in particular on the asym- 
metries present in the flow velocity and density. In an axially-symmetric waveguide, the 
superfluid dissipates through the nucleation of circular ring vortices that shrink across the 
flow within the barrier regions. As soon as the axially symmetry is broken, the rings are 
strongly deformed and can break up into vortex lines or not form at all. The phase slip 
instability observed here can be considered as a superfluid flow dissipation mechanism, since 
it subtracts energy from the potential flow of the superfluid. However, GPE does not in- 
clude the dissipative processes which bring the system to thermal equilibrium. In a realistic 
situation, the vortices, carrying the energy subtracted from the superfluid, are supposed to 
eventually thermalize, with the result of heating the system. The effect of dissipation on 
vortex nucleation has been studied by adding a dissipation term to the GPE 15] . The role 
of finite temperature on the superfluid instability scenario presented here deserves further 
study. 
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